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Aging effects in the two-time correlation function and the response function after a quench 
from a high temperature to some low temperature are considered for a simple kinetic ran- 
dom energy model exhibiting stretched exponential relaxation. Because the system reaches 
thermal equilibrium for long times after the quench, all aging effect are of a transient na- 
ture. In particular, the violations of the fluctuation- dissipation theorem are considered and 
it is found that the relation between the response and the two-time correlation function 
depends on another function, the so-called asymmetry. This asymmetry vanishes in equi- 
librium but cannot be neglected in the aging regime. It is found that plots of the integrated 
response versus the correlation function are not applicable to quantify the violations of the 
fluctuation- dissipation theorem in this particular model. This fact has its origin in the 
absence of a scaling form of the correlation. 

PACS numbers: 05.40.-a, 64.70.Pf, 61.20.Lc 



I. Introduction 

In the last decade, much attention has been paid to the out-of-equlibrium dynamics of 
glassy systems. In particular, the deviations from the fluctuation-dissipation theorem 
(FDT) have been studied for a variety of models[l]. Usually, the extent of the violations 
of the FDT are quantified via the so-called fluctuation-dissipation ratio (FDR) X(t,t w ), 
defined by 

r>U j. \ X(t,t w ) dC(t,t w ) 

R{t f t u ) = — . (1) 
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Here, R(t, t w ) denotes the linear response of a dynamical variable M to a field that has been 
applied to the system at t w and C(t,t w ) is the corresponding auto-correlation function. In 
equilibrium, X(t,t w ) = 1 and the FDT is recovered. In non-equilibrium situations, it is 
tempting to use the FDR for the definition of an effective temperature characterizing the 
state of the systemj2]. For some models of glassy relaxation, it has been found that in the 
scaling regime X(t,t w ) is a function of the correlation alone, X(C). Furthermore, often 
X(C) = 1 holds for short times, whereas X(C) < 1 in the long-time sector, implying an 
effective temperature which is higher than the bath temperature. In addition to a number 
of analytical calculations, several molecular dynamics simulations have been performed on 
model glassforming liquids, for a recent review see ref. [Hj. Of course, in order to allow 
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a meaningful definition of a thermodynamic quantity like a temperature, X(t,t w ) should 
be independent of the dynamical variable used for its calculation. In case that so-called 
neutral variables M[4] are used, X(t,t w ) often is found to be independent of M and in 
some cases the long time limit X°° allows the definition of an effective temperature. 

In a previous paper 0, I considered in detail the relation between the response and the 
correlation for stochastic Markov processes in which case the dynamics is described by a 
master equation A treatment of FDT violations for stochastic models with Langevin 
dynamics has been given earlier [Zj. In both cases it is found that the response is not 
determined by time derivatives of the correlation alone but that an additional function 
A(t,t w ), called asymmetry, is needed to determine the response. In the context of master 
equations, the asymmetry has been discussed for the first time in refs. jSJ E] in the context 
of the aging dynamics in spin glasses. It was shown in ref. that different kinds of 
dynamical variables can give rise to qualitatively different results regarding the asymmetry 
and therefore also the FDR X(t,t w ). Two special choices for the dynamical variable M(t) 
have been considered. One class of variables randomizes completely with any transition 
among the states of the system. Another class of variables is chosen in such a way that 
there is no correlation to the transitions among the states at all. The first type of variable 
is standard in the investigation of trap models fTTH [TTJ EE21 EH]- Even though in both cases 
the variables can be chosen as neutral, the behavior of the relevant dynamic quantities 
can be different. In general, the asymmetry vanishes for the first class of variables under 
very mild conditions but not for uncorrelated variables. In the trap model both kinds 
of variables yield identical results. Furthermore, all dynamical quantities show a (t/t w )- 
scaling in the long-time limit and the FDR X(t, t w ) as determined by eq.(0) coincides 
with the slope in a fluctuation-dissipation (FD) plot of the integrated response versus the 
correlation function |12j. 

In the present paper, I will consider another class of models in which equilibrium is 
reached after long times. Therefore, all aging effects are of a transient nature but the FDT 
still is violated for short enough waiting times. It will be shown that in a simple model 
the asymmetry does not vanish and has a strong impact on the FDR in the aging regime. 
Additionally, it turns out that the FDR calculated directly from its definition deviates 
from the slope in FD plots. After a short review of the calculation of the response and 
the correlation for Markov processes I will discuss the random energy model^l] with a 
kinetic rule that has been considered earlier by Koper and Hilhorst for the relaxation 
in equilibrium. 

II. FDT violations for stationary Markov processes 

In this chapter I briefly recall the derivation of the fluctuation-dissipation relation for 
stationary Markov processes as it has been derived in ref. 5J. I will concentrate on stationary 
Markov processes for simplicity, but also nonstationary processes can be treated in the same 
way. 
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In order to calculate the fluctuation-dissipation relation both, the response and the cor- 
relation are needed. In all following calculations it is assumed that the system is quenched 
from a high temperature to a low 'working' temperature in the beginning of the experimen- 
tal protocol. After a waiting time t w has evolved after the quench, either a field H5(t — t w ) 
is applied and the response is measured or the auto-correlation function is monitored. For 
a dynamical variable M(t) the correlation function is given by: 

C(t,t w ) = (M(t)M(t w )) =J2M k M l G kl (t-t w ) Pl (t w ) (2) 

k,l 

Here, G k i(t — t w ) is the conditional probability to find the system in state 'k' at time 
t, provided it was in state T at time t w . The population pi(t w ) = J2 n Gin(t w )Pn gives 
the probability to find the system in state '1' at the beginning of the measurement. The 
initial high-temperature state of the system is assumed to be described by a fixed set of 
populations, pi — p k (t — 0) with J2 k P k = 1- Furthermore, M k denotes the value of M 
in state 'k'. In the above expressions, I used a discrete notation and the term 'state' is 
meant to represent the stochastic variable under consideration. In many models for glassy 
relaxation, this variable will be a (free) energy. The conditional probability G k i(t — t w ) 
obeys a master equation 6 : 

§;Gki{t) = -£ W nk G M (t) + ]T W kn G nl (t) (3) 

with W k i denoting the rates for a transition from state I to state k. 

For a calculation of the linear response to a field H (t) = H5(t — t w ) conjugate to M, 



R(t, t w ) — 



(4) 

H=0 



5H(t 

one has to fix the dependence of the transition rates on the field H. There is no general 
way to do this and here, as in ref. [S], I use the following form of multiplicatively perturbed 
transition probabilities 

Wif ) = w^eW^M,) (5) 

Here, 7 and fi are arbitrary parameters. If fi + 7 = 1 holds additionally, then the rates 
Wm obey detailed balance also in the presence of the field. 

Without going into the details of the calculation of the response according to eq.(JU), 
I only mention that time-dependent perturbation theory is used to calculate G kl (t) in 
linear order with respect to the field [TS]. In case of stationary Markov processes described 
by a continuous time master equation the result can be cast into the form[S]: 



R(t,t w )=p 



■ dC{t,t w ) _ dC{t,t w ) _ 

7 dt w dt ^ A ^> 



(6) 



with the asymmetry 



A(t, t w ) = M k M t G kn (t - t w ) [Winp n (t w ) - W nm (t w )} (7) 

k,l,n 
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It is important to point out that the asymmetry A(t,t w ) cannot be related to a time 
derivative of the correlation function and therefore the response is not determined by 
C(t,t w ) alone in the general case. If the system is in equilibrium, one has A eg (t,t w ) =0. 
This can be achieved in two ways. Either the system has been prepared in an equilibrium 
state initially, i.e. p° k =p e k q ) or the system has reached equilibrium after a sufficiently long 
waiting time. In both cases one has pi(t w ) =Pi 9 and use of the detailed balance condition 
Wi n p e £ = W n ip e x q shows that A eq (t,t w ) = 0. In these situations C(t,t w ) and R(t,t w ) are 
time-translational invariant and one finds 

R eq (t)=-^ + »)^- (8) 

which for /i = l — 7 is just the well known FDT. 

In the following, I will consider the same types of dynamical variables that have been 
discussed in detail in ref.jHj. One class has been denoted as 'uncorrelated' and the other one 
as 'randomizing'. The meaning of these terms is the following. An uncorrelated variable 
is completely decoupled from the transitions among the states of the system. This is the 
choice made by Koper and Hilhorst^Hl m their treatment of a kinetic random energy 
model. In contrast, a variable for which every transition among the states of the system 
gives rise to a complete randomization is called a randomizing variable. This is the type of 
variable that is used in most calculations for trap models. In the following, I will restrict 
the discussion to variables the values of which are independent of the state of the system, 
M k = M V7c, and furthermore assume distributions of zero mean and unit variance, 

(M) = and (M 2 ) = 1 (9) 

With this choice one explicitly finds for the two types of variables [5 

uncorrelated variables: C(t, t w ) = V] G kk {t — t w )pk{t w ) 

k 

A(t, t w ) = ]T G kl (t - t w ) [W klPl (t w ) - W lkPk (t w )] (10) 
k,i 

and 

randomizing variables: U(t, t w ) = ^ e~ Kk ^~ tw ^p k (t w ) with K k = ^Wi k 

k I 

A(t,t w )=0 (11) 

In the latter case, I denoted the correlation function by U(t,t w ), because this function 
is identical to the probability that the system has not left the state occupied at t w 
considered in the classical treatment of trap models [HH ITTj . For trap models one has 
G kk (t — t w ) = e~ Kk ^ t ~ tw ' ) + 0{N~ 1 ) and therefore uncorrelated variables yield the same 
result as randomizing variables for large N. 

Before closing this section, some comments are in order. One comment concerns the 
values of the parameters 7 and /1. Even if one enforces detailed balance to hold in the 



4 



presence of the field, i.e. if 7 + fi = 1 holds, there is no general way to fix them. Only 
in some special cases one has some arguments based on the physics of the corresponding 
model. For instance, if one considers particles hopping on a (disordered) lattice, one would 
naturally choose 7 = /x = 1/2. Slightly more generally, one would expect the choice 
7 = /i = 1/2 to be meaningful if the master equation considered has a form allowing a 
Kramers-Moyal expansion^. These arguments, however, are not sufficient to allow the 
determination of 7 and /x in general. 

Regarding the asymmetry A(t,t w ), it was stated above that this quantity cannot be 
related to a time derivative of the correlation function. Furthermore, there is no obvious 
physical interpretation of A(t,t w ). Apparently, it has to be related to some correlations 
among the dynamical variables as otherwise it would not vanish in case of randomizing 
variables. However, A(t,t w ) appears quite strongly to depend on both, the model consid- 
ered and the choice of variables. At present, it is not clear whether it is related to some 
other dynamical function of general importance. 

III. A kinetic random energy model 

In the following, I will apply the results of the preceeding section to a very simple model 
which allows to calculate all quantities analytically. The model exhibits interrupted aging 
only, reaching equilibrium in the long-time limit. Such kind of models are interesting to 
study because in canonical glasses the situation is somewhat similar. If a glass-forming 
liquid is quenched to a temperature below the glass transition temperature, this means that 
the re-equilibration time exceeds typical experimental time scales. However, one expects 
the system to reach (metastable) equilibrium for very long times. A similar situation may 
appear in computer simulation studies of model glass-formers. Here, the available computer 
time sets the time scale for following the equilibrium properties of the system. 

I will treat a special case of a kinetic random energy model^l], following the treatment 
of Koper and Hilhorst|15j. who considered various forms of the transition rates. For the 
purpose of the present discussion the simplest choice 



is sufficient, where kq is a rate constant, to be set to unity in the following, /3 = T _1 , and 
B k is a Boltzmann factor. According to eq. (fT2"|) the W k i only depend on the destination 
state of the transition and the corresponding ME can easily be solved with the result 



with the partition function Z(/3) = J2k Bk- The equilibrium populations are given by 
p1 q = Gik{t—>- 00) = Z(f3)~ 1 Bi. Therefore, the system reaches equilibrium for long times and 
all aging phenomena are of a transient nature. 

As already noted, I will use a distribution of 'magnetizations' M with (M) = and 
(M 2 ) = 1 throughout the calculations. In addition, the initial populations p° will be 



Wm = K-oB k = K e 



(12) 



G lk (t) = Z(py 1 B l + 5 lk -Z(P)~ l B l exp(-Z{p)t) 



(13) 
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chosen as pf = iV _1 , appropriate for a quench from high temperatures in the beginning 
of the experimental protocol. The final temperature is chosen to be lower than the phase 
transition temperature T c of the random energy model Jl] . The calculation of the quantities 
of interest proceeds in the following way. First one calculates the correlation and asymmetry 
according to eqns. ljlOjllj) using the expression given above for the Greens function Gkiif) 
and Pkitw) = J2iGki{t w )p*i. The quantities calculated this way depend on the Boltzmann 
factors, which will be denoted by a superscript B, e.g. C^ B \t,t w ) etc.. These expressions 
have to be averaged over the distribution of Boltzmann factors, which derive from the 
distribution of the random energies. Below T c the random energies are exponentially 
distributed and therefore the corresponding distribution of Boltzmann factors is given 
by p(B) = (v/N) x B- X ~ x for (v/N)^ < B < oo and p(B) = otherwise JEj. Here, 
x = T/T c , v = (2a/ vr log 2)" 1 and iV denotes the number of random energies, to be sent 
to infinity at the end of a calculation. In a shorthand notation one then has F(t,t w ) = 
(F^(t,t w )} = JdBp(B)F^(t,t w ). In the following, I will discuss the two classes of 
variables introduced above separately. 



Uncorrelated variables 

For this case one finds, using the abbreviation Z{(3) = Z{2(3) / Z{(3) 2 : 

C( B \t, O = Z(P) + [1 - Z{(3)\ exp (-Z(P)(t - t w )) 

+Z{p) [exp (-Z(J3)t) - exp {-Z{(5)t w )\ (14) 

and 

A^ B \t,t w ) = Z{f3)Z{f3) [exp (-Z((3)t w ) - exp (-Z(f3)t)} (15) 

This expression can be cast into the form A^(t,t w ) = d t C {B \t,t w ) + d tw C {B \t,t w ), which 
allows to write the expression for the response in the form: 

i**>(t,O = -/3(7 + A0 dCiB) g f tw) (16) 

Therefore, in this case the response depends on the parameters 7 and \x only via their 
sum. As already noted, one has 7 + fi = 1 if it is assumed that detailed balance holds 
in the presence of the field. In order to perform the averages of the quantities given in 
eqns. (jl4ll5H6|) . the only integrals required are (e~ z ^ t ) and (Z(P)e~ z ^ t ) . The calculation 
of these integrals has been performed in ref.p3] with the result: 

$(t) = ( e - z{li)t ) = exp (-vt x ) 

#(t) = (Z{(3)e- z{/3)t ) = (1 - x) [$(t) - v l/x tT{\ - l/x; vt x )] (17) 

where v = vT(l — x) with the Gamma function T(a). Furthermore, T(a;b) denotes the 
incomplete Gamma function. In terms of these functions the averaged quantities read: 

C(t, t v ,) = l-x + $(t- t w ) -y(t-t w ) + - V(t w ) 
A(t,t w )=d t *(t)-d Ue *(t w ) (18) 
R(t,t w )=-(3 (j + n)d t C(t,t w ) 
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As stated above, the system reaches equilibrium for long waiting times t w 3> r eq , where r eq 
denotes the re-equilibration time. In this case, one finds from eq. (jl8|) . assuming r = (t — t w ) 
to be finite: 

Ceq ( T ) = 1 - x + $(r) - *(r) and it! eq (r) = -/? (7 + /i) (19) 

and Aeq(t,t TO ) = 0, as expected. Note that C(t,t w ) can be cast into the form C(t,t w ) = 
C eq (t — t w ) + C ag (t,t w ). However, C ag (i, i w ) does not exhibit a (t/t w ) -scaling. 

The behavior of C eq (r) has been discussed in detail by Koper and Hilhorst [T3] . C eq (r) 
decays in a stretched exponential fashion, C eq (r)=C^ + (1 — exp { — (£/r eq )^ cq } with 
Cg? = (1— x). Also for finite t w C(t,t w ) can be fitted to a Kohlrausch function C(t w +r, t w ) = 
C°°(t w ) + (1 — C°°(t w )) exp { — (t/t k ) i3k } where the plateau-value C°°(t w ) equals zero for 
vanishing t w and smoothly approaches for long waiting times. This behavior is shown 
in Fig. la, where C(t w + r,t w ) is plotted versus r for various waiting times (upper panel). 
The dotted lines represent least-square fits to the mentioned Kohlrausch function. Only for 
intermediate t w there are some minor deviations visible. In the lower panel of Fig. la the 
correlation function is shown in a normalized way so that it decays from unity to zero. It is 
evident from that plot that the relaxation time for intermediate waiting times exceeds the 
equilibrium relaxation time. This fact is further substantiated in Fig. lb, where the results 
of Kohlrausch fits to C(t w + r, t w ) are shown. The relaxation time TK{t w ) first increases and 
then decreases again before it tends to the limit r eq . The stretching parameter also first 
increases and goes through a maximum before it reaches its equilibrium value. This latter 
fact means that the distribution of populations first narrows and then broadens again as a 
function of the waiting time. A qualitatively similar behavior has also been found in some 
computer simulations |18j and in a free-energy landscape model for glassy relaxation |19j. 
Another important feature that becomes evident from Fig.l is that the correlation function 
does not exhibit a -scaling, cf. the lower panel of Fig. la. For other temperatures the 
behavior is the same. 

Next, the asymmetry is considered. According to eq.([18|). A(t,t w ) reaches a plateau 
determined by [— dt w ^(t w )] for long times t and finite t w . Furthermore it can be shown, 
that A(t w + r, t w ) ~ r for small r. The asymmetry A(t w + r, t w ) is plotted versus t in Fig.2 
for various values of t w for x = 0.3. The behavior for other values of x is very similar. The 
linear behavior at short times is seen in the lower panel, where a logarithmic scale is used. It 
is evident that the cross-over from the linear behavior to the plateau takes place around t w , 
t ~ t w . It is obvious from Fig.2 that for uncorrelated variables the asymmetry significantly 
contributes to the response for short waiting times and thus cannot be neglected. 

According to eqns.Q and ()10|) the FDR in the present case is given by: 

L7 + /XJ A (t, t w ) - dtwCM + mt _ tw) _ mt _ tw) + Qt ^ {tw) W 
This FDR may be considered as a function of either r = (t — t w ) or as a function of t w . In 
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the first case, one finds the following limits (/x = l — 7): 

X(t w + r, t w ) -> 1 ; r < t w and X(t w + r, t w ) -> ; r > t w (21) 

Similarly, if X(t,t w ) is considered as a function of t w , as has been suggested in ref. [I], one 
finds 

X(t,t w ) ^1 ; t w ^t and X(t, t w ) -> ; t w -> (22) 

Furthermore, one finds from the limiting behavior of the function ^(t) that X(t, t w )~t]^ x 
for small t^,. 

This behavior of the FDR is illustrated in Fig. 3 for 2 = 0.3. In the upper panel X(t w + 
T,t w ) is plotted versus r for several values of the waiting time t w . It is seen that the 
crossover from X — 1 to X = takes place around r~t w . This means that for finite t w one 
has X°°(t w ) = 0. The lower panel of Fig. 3 shows X(t,t w ) as a function of t w for several 
values of the later time t. It is evident how X approaches the equilibrium value X = 1 for 
t w —*t. The dotted line shows the mentioned t^ _:r -behavior for t w <^t. 

For the integrated response x{t,t w ) = J t * dtfR(t,t') one finds: 

X (t, t w ) = /% + fi)[x- $(t - t w ) + V(t - t w ) -(t- t w )d t ^(t)] (23) 

where the first three terms correspond to the equilibrium response, Xeqi^) — (3 [1 — C e0L {r)] 
(because of C eq (0) = 1), which is reached for t w ^> r eq . The limiting value for long times 
is thus given by Xeq(oo) =(3x = T~ . For finite t w , the limiting behavior of the integrated 
response is given by xi^w + t, t w ) — > (7 + ^T^ 1 for r — > 00, which for /x = 1 — 7 coincides 
with the equilibrium value. 

In Fig. 4, typical FD-plots, /3 _1 x(^o + T ^w) versus C(t w + r,t w ), for x = 0.3 (upper 
panel) and x = 0.6 (lower panel) are shown. Here, r is used as the curve parameter and 
the different curves are for various t w - values. Note that C(t,t) = 1 always holds in this 
model. Furthermore, it is assumed that fi = 1 — 7 (otherwise there would be an extra factor 
7 + /i). The dotted lines represent slopes of (—1) (FDT) and of (—x). From the FD-plots, 
for t w < r eq (r eq (x = 0.3) ~ 80 and r eq (x = 0.6) ~ 2.5), one would naively extract a limiting 
FDR X°° that coincides with the transition temperature T c . For longer waiting waiting 
times, the initial slope of (—1) persists for longer times and the FDT is recovered. 

From the discussion of the FD-plot and the FDR it is obvious that the slope in a FD- 
plot is not related to the FDR in this simple model. This fact may be of importance if 
one cannot perform experiments over a long enough time scale in order to monitor the 
re-equilibration of the system. Note that one might argue that the FD-plots built from a 
correlation function with a finite long-time plateau value may give rise to erroneous results. 
However, in the present case one has C°°(t w ) ~0 for small t w , cf. Fig. la. 

Randomizing variables 

For this choice of variables one finds from eqn. (jllj) after averaging over the random energies: 

Il(t,t w )=®(t-t w ) (24) 
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where $(£) is given in eq.(|17|). This means that the correlation and therefore also the 
response are time-translational invariant and for 7 + fi = 1 the FDT holds because of 
A Tand (t,t w ) = 0, cf. eq. lfTTj) . Thus, randomizing variables do not at all allow to probe 
the out-of-equilibrium situation imposed by the quench from high temperatures into the 
low temperature phase for this simple model. This can be understood from the fact that 
the dynamics is not thermally activated for the simple choice (112)1 of the transition rates. 
Therefore, the variable has lost memory of the initial non-equilibrium situation completely 
after a single transition. This is in contrast to the situation in the trap model, where 
aging is due to the fact that a equilibrium distribution of populations does not exist in the 
low-temperature phase [T7|. 

V. Conclusions 

In the present paper I have recalled the results obtained earlier for the relation between the 
correlation function and the linear response of a system exhibiting stochastic dynamics de- 
scribed by a stationary continuous time Markov process. For the calculation of the response 
the transition rates of the master equation are assumed to be perturbed multiplicatively 
by the applied field. In general, the response is not determined by time-derivatives of the 
correlation alone, but the asymmetry A(t,t w ) has to be considered additionally, which in 
general cannot be related to time derivatives of the correlation function. This function in 
general only vanishes under equilibrium conditions, in which case the FDT is recovered, if 
the dependence on the field of the transition rates is chosen in a symmetric way (7 + /i = l). 

I considered two different classes of variables, namely variables that are uncorrelated 
with the states of the system and variables which randomize completely with each transi- 
tions among the states. Both types of variables are neutral in the sense that their values 
are not correlated with the state of the system, in particular the energy 

As an example for a simple model of glassy relaxation I considered a kinetic REM with 
an extremely simple choice for the transition rates. This model shows intermittant aging 
behavior and reaches equilibrium for long times. If uncorrelated variables are chosen, 
the asymmetry is finite for short waiting times. Furthermore, in this case the limiting 
slope in an FD-plot is determined by the transition temperature of the model. However, 
this slope does not coincide with the value of the FDR as determined from its definition. 
Therefore, in this case a FD-plot does not yield the correct value for the FDR. Of course, 
the reason for this finding lies in the fact that the correlation does not obey a scaling-law 
but rather decays in a stretched exponential manner. Even more relevant to experimental 
determinations of the FDR in glassy systems may be the fact that the results for different 
variables are different. If randomizing variables are used instead of uncorrelated variables, 
the asymmetry vanishes and all quantities are time-translational invariant. Thus, one 
always has a FDR that equals unity in this case. 

To conclude, I have shown that for some general class of models the fluctuation- dissipa- 
tion relation is determined by time-derivatives of the correlation function and an additional 
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function, the asymmetry A(t,t w ). In general, the asymmetry can be shown to vanish for 
randomizing variables under the conditions considered usually. For the model considered 
in this paper, A(t,t w ) has a strong impact on the behavior of the FDR for uncorrelated 
variables. Additionally, it appears that one must be careful in the determination of the 
FDR X(t,t w ) from FD-plots in some cases. 
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Figure captions 



Fig.l : a) The correlation function for an uncorrected variable in the kinetic REM 
plotted as a function of time for x = T/T c = 0.3 and various values of the waiting time 
t w . Upper panel: C(t w + r,t w ) for \og 10 (t w ) = —5,-2,0,1,2,5. The dotted lines 
represent fits to a Kohlrausch function. Lower panel: C(t w + T,t w ) normalized in a 
way that it decays from unity to zero. This plot shows that for intermediate waiting 
times the relaxation time goes through a maximum. 

b) Fitting parameters t k (tw) (upper panel) and /3^(t w ) (lower panel) resulting from 
Kohlrausch-fits C(t w + T,t w ) = C°°(t w ) + (1 - C°°(t w )) exp {-(r / W K f K } ■ 

Fig. 2 : A(t w + r, t w ) versus r, for x = 0.3 for various waiting times t w for the kinetic REM 
with an uncorrelated variable. Upper panel: linear scale, lower panel: logarithmic 
scale. In the lower panel, the dotted line is proportional to r. 

Fig. 3 : Upper panel: X(t w + T,t w ) versus r for x = 0.3 and \og 10 (t w ) = —5,-2,0,2,5 
for the kinetic REM with an uncorrelated variable. The dotted vertical lines show 
the corresponding values of t w . Lower panel: X(t,t w ) versus t w for x = 0.3 and 
log 10 (t) = —5, —3, —1, 1, 3, 4, 5. The dotted line is oc t]^ x . 

Fig. 4 : fi^xitw + T,t w ) versus C(t w + r, t w ), for x = 0.3 (upper panel) and x = 0.6 
(lower panel) for various waiting times t w for the kinetic REM with an uncorrelated 
variable. The waiting times chosen are \og w (t w ) = —5,-1,0,1,2 for x = 0.3 and 
\og 10 (t w ) = —5,-1,0,2, for x = 0.6. The dotted lines represent the slopes expected 
for the FDT (slope: —1) and a slope of —x. 
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G. Diezemann; Fig. 1b 
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G. Diezemann; Fig. 3 
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